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ABSTRACT 


Presently, the United States Navy is searching for an 
improved method to predict the damage to a ship hull or 
underwater structure that results from an underwater 
explosion. One method of predicting this damage is through 
the use of nonlinear finite and boundary element analysis. 
Underwater Shock Analysis (USA) code combined with 
VEC/DYNA3D code is used for the analysis of the effect of 
explosive shock on numerical models. Initial geometric 
[Mileerrections are introduced in the numerical model using 
modal imperfections. The resulting numerical model is then 
subjected to a simulated underwater shock using the combined 
USA/DYNA3D code. A sensitivity analysis is performed to 
look into the details on the damage resulting from these 


Smmhations. 
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I. INTRODUCTION 


The use of USA/DYNA3D for analyzing structures 
Subjected to underwater shock has been shown to be an 
effective tool in predicting the response of these 
Structures and the resulting damage (Chisum 1992). In some 
cases the actual deformation of a cylinder subjected to an 
underwater shock 1S not accurately predicted. One 
Seeleomatlom for “this inaccuracy 1s that actual cylinders 
have many imperfections (e.g., out of roundness, thin 
sections, voids, etc.) and the finite element modeling of 
m@ege Cylinders often do not take into account the initial 
imperfections ice are present. By THe rocucing 
imperfections in the position of the node points generated 
by a finite element mesh generator it is possible to more 
accurately model the geometry of an actual cylinder. With a 
more accurate model of the actual cylinder the finite 
element numerical analysis of the resulting damage due to 


underwater shock may be more accurate as well. 


II. INITIAL IMPERFECTIONS 


The location of the node points of a structure modeled 
by a finite element mesh generator are precisely located (to 
the numerical accuracy of the modeling code) at the 
positions specified in the inputs to the mesh generator 
alee ie min. Manufactured structures usually have many 
imperfections that are not accurately modeled by many finite 
element mesh generators. Surveys have been performed to 
measure the imperfections that Natura ly. exist aia) 
cylindrical shell structures (Arbocz 1982). These data have 
been collected into data banks, and Arbocz suggests that 
they be used to improve design criteria for buckling of than 
shells. The imperfection data banks show that imperfections 
have characteristic distributions that include decreasing 
modal amplitudes with increasing mode number. Kirkpatrick 
(1989) found that by introducing initial modal imperfeceweme 
in a cylindrical structure subjected to blast shock loadaimg 
that the resulting numerical analysis agreed much closer to 
the experimental data than the analysis without these 
initial imperfections. 

Initial imperfections were introduced into the location 


of the node points of the finite element model using a 


Summation of modal imperfections expressed as the cosine 


series shown in Eq. (1). 


N 
AR(8) = >) Ancos(nO + on) pe ee tyes 20> ct a Bi 


n=2 


where AR is the radial imperfection, 9 is the angular 
position, N is the maximum modal contribution, An is the 
modal amplitude, and ¢ is a random modal phase shift. The 
assumption that the modal phase shift is a random variable 
1s reasonable for many shells. An empirical form for the 


modal amplitude is expressed as shown in Eq. (2). 





where An 18S the modal amplitude of the ath modal 
imperfection, nm 1s the mode number and xX ander _ are 
coefficients used to fit the available data for shells of a 
given construction. The modal amplitude can be modeled as a 
constant percentage of the shell prmremnecas For many of the 
Studies in this investigation the modal amplitude is 
Seleeant as shown in Eq. (3) at 1% of the shell thickness 


feneall modes. 


n=Q.01h i Ope wers ILM 7: Ce gees ee a a) 


Kirkpatrick found that the modal amplitudes in Eq. (4) and 


(5) accurately modeled the imperfections of his test 
cy lander: 
An = 0:05/: NSO on ee ee 
2h 
An=— n27 (59) 





The modal imperfections introduced in this study are 
assumed to be a function only of the circumferential 
position of the node and not a function of the axial 
position of the node. The programs used to modify the node 


positions are shown in Appendix A. 


III. TWO DIMENSIONAL INFINITE CYLINDER MODELS 


FOX (1992) developed a two dimensional infinite 
cylinder model to validate the use of the USA/DYNA3D code 
for use in underwater shock analysis. This model was 
further developed into the full symmetry two dimensional 
infinite cylinder model shown in Figure 1. This model has 
40 elements around the circumference of the cylinder, radius 
of 6 inches, shell thickness of 0.06 inches and an element 
meaegiem OF 0.006 inches. The same restriction on the length 
of the elements found for the model Fox developed 
(length/radius=0.001) was also true for this model. Ox 
Stated that longer elements resulted in "residual three 
dimensional effects" that caused oscillations in the 
numerical results. Symmetric boundary conditions were 
imposed on both ends of the cylinder to model an infinite 
length cylinder. This model and all subsequent models 
(unless specifically noted) use Belytschko-Tsay (Belytschko 
1984) shell elements for their numerical efficiency and to 
eliminate the problems of using a solid element to model a 
thin walled structure. For details of the model input 
parameters see Appendix B. The shell material is modeled as 
mild steel, an elastic-plastic material with a Young's 


feats sof 2,010’ psi, a Poisson's ratio of 0.3, a 


Figure 1. Two Dimensional Infinite Cylinder Model 


hardening modulus of 5100 psi anda yield stress of 32,000 
psi. For details of the material models used in this study 


see Appendix C. 


A. Perfect Cylinder Model 

This cylinder was subjected to a underwater shock wave 
modeled as a plane wave with an amplitude of 500 psi anda 
duration of 1 millisecond resulting from an explosive charge 
located along the x-axis. The pressure time history of the 
shock wave is shown in Figure 2. The early damage pattern 
of the cylinder at 0.5 milliseconds is shown in Figure 3. 
The front face of the cylinder is flattened by the impacting 
plane wave approaching from the right. However, Figure 4 
shows that the final damage pattern at 5 milliseconds is 
very different than this early deformation pattern. Pak 
plastic deformation has been completed by 5 milliseconds. 
The side of the cylinder facing the approaching shock wave 
shows a pronounced protrusion. This pattern was unexpected 
based on experimental results which have usually shown 
inward deformations of a cylindrical shell facing the 
explosive charge. It 1s hypothesized that this raised 
section 1s due to the fact that the modeled cylinder is a 
perfect cylinder and is not representative of a true 


manufactured cylinder. The manufactured cylinder will have 
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0.5 ms (deformations scaled up 100X) 


Figure 3. Deformation of Two Dimensional Infinite Cylinder 
Subjected to a Square Pressure Pulse of 500 psi for 1 ms 
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5.0 ms (deformations scaled up 10X) 


Figure 4. Deformation of Two Dimensional Infinite Cylinder 
Subjected to a Square Pressure Pulse of 500 psi for 1 ms 


many imperfections in the shell which is not present in the 
computer models. 

In order to study the effect of mesh size on the 
response of this cylinder another model was developed with 
64 elements around the circumference of the cylinder. 
Figure 5 shows that the response of this cylinder is similar 
to that for the 40 element Cy lange. There is still a 
distinct protrusion of the shell toward the direction of the 
approaching shock wave. All subsequent cylinders will use 


40 elements around the circumference of the cylinder. 


B. Imperfect Cylinder Models 

The introduction of initial imperfections signifive@amme, 
changes both the shape and magnitude of the resulting 
deformation. Introduction of a 6th mode imperfection, as 
shown in Figure 6 (imperfection magnitude scaled by a factor 
of 100) with a modal amplitude of 5% of the shell thickness 
and no random phase shift, results in the damage pattern 
Shown (deformations scaled by factor of 10) in Figure 7. A 
mode 6 imperfection was chosen because of the ease in 
identifying. 26s) distime eivew saetcmnr Use of another model 
imperfection would only have changed the shape of the final 
deformation not any conclusions resulting from this study. 
Comparison of the final damage pattern with the shape of the 


initial imperfection shows that the damage pattern follows 
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5.0 ms (deformations scaled up 10X) 
Figure 5. Deformation of Two Dimensional Infinite Cylinder 


Subjected to a Square Pressure Pulse of 500 psi for 1 ms 
Fine Mesh 


1] 


Initial Mode 6 Imperfection (imperfections scaled up 100X) 
Figure 6. Two Dimensional Infinite Cylinder Subjected to a 


Square Pressure Pulse of 500 psi for 1 ms, Mode 6 Initial 
Imperfection As=0.05h no phase shift 


Shock 
Wave 


5.0ms (deformations scaled up 10X) 
Figure 7. Two Dimensional Infinite Cylinder Subjected to a 


Square Pressure Pulse of 500 psi for 1 ms, Mode 6 Initial 
Imperfection As=0.05h no phase shift 


2 


the initial imperfection resulting ina final damage pattern 
looking like the mode 6 initial imperfection magnified due 
to the external shock wave applied to the cylinder. The 
introduction of imperfections caused the deformation of the 
cylindrical shell to follow the initial imperfection pattern 
of the shell. In addition to changing the damage pattern, 
the initial imperfection also causes much greater 
deflections of the shell for the same shock wave intensity. 
Figure 8 shows the deformation of a two dimensional cylinder 
with an initial imperfection 5% of the shell thickness with 
emeeea mode 10 contribution and no random phase shift. 
Comparison of Figures 7 and 8 show that increasing the mode 
number of the initial imperfection while keeping the 
amplitude constant at 5% of the shell thickness results in 
greater damage to the cylinder for the same shock wave 
pressure. This is due to the greater local curvature of the 
shell with imperfections of higher mode number than for 
lower mode numbers. Figure 9 shows how the total internal 
strain energy (a measure of the damage to the cylinder) of 
the cylinders increases as the mode number of the initial 
imperfection increases while the amplitude of the 
imperfection is constant at 5% of the shell thickness for 
the same shock wave pressure. 

Another cylinder was modeled with aimperfections 


containing the first 10 modal imperfections, as shown in 


ae 


4 


5.0 ms (deformations scaled up 3X) 


Figure 8. Two Dimensional Infinite Cylinder Subjected to a 
Square Pressure Pulse of 500 psi for 1 ms, Mode 10 Initial 
Imperfection A1io=0.05h no phase shift 
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Figure 9. Total Internal Strain Energy for Cylinders 
Subjected to a Square Pressure Pulse of 500 psi for 1 ms 
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Figure 10 with modal amplitudes for each mode of 1% of the 
shell thickness and with random phase shifts. The resulting 
damage pattern is shown in Figure 11. Again, the general 
Shape of the final deformation pattern is very close to the 
shape of the initial imperfection pattern. 

The introduction of a different initial imperfection, 
shown in Figure 12 with initial imperfections magnified by 
50 times, using the modal amplitude of equations (4) and 
Oyeeeresults in the final deformation pattern shown in 
Figure 13. Again, the deformation of the shell of the 


cylinder follows the initial imperfection pattern. 


Initial Imperfection (imperfections scaled up 200X) 
Figure 10. Two Dimensional Infinite Cylinder Subjected to a 


Square Pressure Pulse of 500 psi for 1 ms, Initial 
Imperfection First 10 Modes An»=0.01Ih with Random Phase Shift 
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5.0 ms (deformations scaled up 2X) 


Figure 11. Two Dimensional Infinite Cylinder Subjected to a 
Square Pressure Pulse of 500 psi for 1ms, Initial 
Imperfection First 10 Modes An»=0.01Ih with Random Phase Shift 


att iad Impertection Prasat ee Scaled up 50x) 


peogeee 12. Two Dimensional Infinite Cylinder Subjected to a 
Square Pressure Pulse of 500 psi ims, Initial Imperfection 


First 10 modes An=0.05h for Modes 2 through 6, An=2h/n’ for 
Modes 7 i ae with Random Phase Shift 


<Shock 
Wave 
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5.0 ms (deformations scaled up 2X) 


Figure 13. Two Dimensional Infinite Cylinder Subjected to a 
Square Pressure Pulse of 500 psi 1ms, Initial Imperfection 


First 10 modes An=0.05h for Modes 2 through 6, An=2h/n’ for 
Modes 7 through 10 with Random Phase Shift 
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IV. THREE DIMENSIONAL RING STIFFENED INFINITE CYLINDER 


MODELS 


The three-dimensional ring stiffened infinite cylinder 
model shown in Figure 14 was developed from the two- 
dimensional infinite cylinder model. This model has a total 
of 400 elements with 40 elements in the circumferential 
direction and 10 elements along the length. Again symmetric 
boundary conditions were imposed on both ends of the 
cylinder™ to model-> van antiniee. Bey ingen. This three- 
dimensional model has stiffeners 0.12 inches thick and 1 
inch deep located on 12 inch spacing. The shell is 0.06 
inches thick. For details in the model input parameters see 
Appendix D The shell and stiffeners are modeled as mild 


steel with the same properties as listed earlier 


A. Perfect Cylinder Model 
This cylinder was subjected to the same plane shock 
wave as the previous models. The early deformation of the 
cylinder at 0.2 milliseconds (Figure 15) shows the outer 
shell pinching in on either side of the stiffener located in 
the center of the cylinder. Figure 16 shows the flattening 
of the face of the cylinder facing the explosive charge at 


O22 ima I iteeeconace At this time in the deformation history 
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Three Dimensional Ring Stiffened Infinite 


migure 14. 


Cylinder Model 
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Front View @ 0.2 ms (deformations scaled up 50X) 


Figure 15. Deformation of Three Dimensional Ring Stiffened 


Infinite Length Cylinder Subjected to a Square Pressure 
Pulse of 500 psa ter ims 
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Wave 





End View @ 0.2 ms (deformations scale up 50X) 
Figure 16. Deformation of Three Dimensional Ring Stiffened 


Infinite Length Cylinder Subjected to a Square Pressure 
Pulse of 500 pSatonr se ms 
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there 1s no evidence of the formation of a local raised area 
on the front of the cylinder facing the charge. IHek= ac bagel 
deformation pattern shown in Figure 17 shows pinching of the 
shell on either side of the stiffener. This pinching effect 
was also noted by Chisum (1992) and appears to be present in 
many finite element simulations of underwater shock which 
model test cylinders as perfect cylinders. Again, the final 
deformation pattern of the cylinder side facing the charge 
has a local protrusion of the shell material toward the 
explosive charge as shown in Figure 18. PNG lal gr wes coampalgs: 
hypothesized that this damage pattern is the result of the 
model cylinder being a perfect cylinder while test cylinders 


will always have imperfections in the shell of the cylinder. 


B. Imperfect Cylinder Models 

The introduction of initial imperfections significantly 
changes both the shape and magnitude of the final 
deformations seen in this cylinder model. The 5% mode 6 
imperfection shown in Figure 19 was introduced into this 
model. The resulting deformation due to a 500 psi 1 
millisecond plane wave pressure pulse is shown in Figures 20 
eoughn 23. The initial imperfection shape is clearly 
evident in the final damage pattern. However the shape of 
the cylinder at 0.2 milliseconds does not show any of the 


Hittial imperfection. There is insufficient time elapsed 


vA 





Front View @ 5.0 ms 


Figure 17. Deformation of Three Dimensional Ring Stiffened 
Infinite Length Cylinder Subjected to a Square Pressure 
Pulse of 500 psi for 1 ms 
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End View @ 5.0 ms 


Figure 18. Deformation of Three Dimensional Ring Stiffened 
Infinite Length Cylinder Subjected to a Square Pressure 
Pulse Gf 5000 psa om eens 
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Mode 6 Initial Imperfection As6=0.05h no phase shift 
(imperfections scaled up 100X) 


Figure 19. Deformation of Three Dimensional Ring Stiffened 


Infinite Length Cylinder Subjected to a Square Pressure 
Pulse of 500 psi for 1 ms 
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(deformations scaled up 50X) 


Front View @ 0.2 ms 
Deformation of Three Dimensional Ring Stiffened 


Figure 20. 
Infinite Length Cylinder with a Mode 6 Initial Imperfection 
Subjected to a Square Pressure 


Ab=0.05h no phase shift, 
Pulse vor 1500 ca tow deans 


¢Shock 
Wave 





(deformations scaled up 50X) 


End View @ 0.2 ms 
Deformation of Three Dimensional Ring Stiffened 


Pigure Zi. 
Infinite Length Cylinder with a Mode 6 Initial Imperfection 
As=0.05h no phase shift, Subjected to a Square Pressure 
Pulse of 500 psi for 1 ms 
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Front View @ 5.0 ms 


Figure 22. Deformation of Three Dimensional Ring Stiffened 
Infinite Length Cylinder with a Mode 6 Initial Imperfection 
Aco=0.05h no phase shift, Subjected to a Square Pressure 
Pulse of 500 psi for 1 ms 
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Figure 23. Deformation of Three Dimensional Ring Stiffened 
Infinite Length Cylinder with a Mode 6 Initial Imperfection 
Acb=0.05h no phase shift, Subjected to a Square Pressure 
Pulse of 500 psi for 1 ms 
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for the deformations to grow large enough to show this 
initial imperfection. However the pinching of the shell on 
either side of the stiffener has already begun at 0.2 
Mid tseeonds: The final deformation still shows this 
pinching effect however the magnitude of this type of 
deformation is reduced by the initial imperfection. 

With the introduction of an initial imperfection Cage 
first 10 modes with modal amplitudes of 1% of the shell 
thickness and with random phase shifts as shown in Figure 24 
the resulting damage pattern again shows thai the 
deformation of the shell will preferentially follow the 
initial imperfections. The resulting damage pattern is very 
different from the perfect case (Figures 25 and 26). Most 
noteworthy is the elimination of the pinch in the outer 
Shell on either side of the stiffener. In addition the 
magnitude of the deformation in the outer shell is immeg 
greater than for the perfect cyiimider) “case: The 
introduction of these imperfections results in sites where 
the cylinder preferentially deforms during the shock 
pressure. Thus when the cylinder deforms the cdeformagzea 
follows the initial imperfections resulting in the final 
Shape of the cylinder looking like the initial imperfection 


Shape. 


— 


ee 


Initial Imperfection (imperfections scaled up 200X) 


Figure 24. Three Dimensional Ring Stiffened Infinite Length 
Cylinder Subjected to a Square Pressure Pulse of 500 psi for 
1 ms, Initial Imperfection First 10 Modes A»a=0.01h with 
Random Phase Shift 
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Front Viewac Ss omeims 


Figure 25. Three Dimensional Ring Stiffened Infinite Length 
Cylinder Subjected to a Square Pressure Pulse of 500 psi for 
1 ms, Initial Imperfection First 10 Modes An»=0.01h with 
Random Phase Shift 
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End View @ 5..0 ms 


Figure 26. Three Dimensional Ring Stiffened Infinite Length 
Cylinder Subjected to a Square Pressure Pulse of 500 psi for 
1 ms, Initial Imperfection First 10 Modes An=0.01h with 
Random Phase Shift 
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C. Parametric Study on the Effect of Cylinder Geometry 
A parametric study changing basic cylinder geometry was 
done to see the effect of these changes on the response of a 
three dimensional ring stiffened infinite cylinder to 


underwater shock. The parameters were varied as follows: 


Shell Thickness Gm@ Gro -0 b2 wand 0.247 1nches 
Stiffener Spacing iin; =o.0 “and (6.0 1nehes 
Diameter 8.0, 12.0 and 18.0 inches 


This geometric variation covers the range of shell thickness 
tO radius ratios from 0.0067 to 0.06 and stiffener spacing 
Memeeaadius ratios from 0.67 to 3.0. 

Increasing the shell thickness strengthened the 
cylinders and in order to obtain significant permanent 
deformation of the shell i1t was necessary to increase the 
amplitude of the input shock wave. The amplitude of the 
Shock wave was adjusted to cause damage to the cylinder 
without completely collapsing the cylinder. This <TR peut 
shock wave pressure ranged from 500 psi for 1 ms for the 
thinnest shell to 2400 psi for 1ms for the thickest shell. 

The resulting damage patterns for a perfect cylinder 
8.0 inches in diameter (shell thickness to radius ratios 
from 0.015 to 0.06 and stiffener spacing to radius ratios 
from 1.5 to 3.0) are shown in Figures 27 and 28. The shock 
wave pressure for the 0.06 inch shell thickness is 600 psi 


for 1ms, for the 0.12 inch shell thickness the pressure 1s 
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1200 psi for 1 ms and for the 0.24 inch shell thicknesome— 
pressure 1S 2400 psi for 1 ms. 

The front view shows that increasing the shell 
thickness causes the pinching of the shell near the 
stiffener to decrease but it is not entirely eliminated at a 
shell thackness of 0.24 inches. Decreasing the stiffener 
Spacing to 8.0 inches from 12.0 inches also has little 
effect on the pinching of the shell but at a stiffener 
Spacing of 6.0 inches there is very little pinching of the 
shell near the stiffener. 

The end view shows that at a shell thickness of 0.06 
inches there is some damage to the front and back of the 
cylinder but most of the damage occurs to the top and bottom 
of the cylinder. For a stiffener spacing of 12.0 inches and 
a shell thickness of 0.12 inches there is still sign2 fieeme 
damage to the cylinder at the top and bottom but there is 
additional damage to the front and back of the cylinder. 
Finally at a stiffener spacing of 12.0 inches and a shell 
thickness of 0.24 inches most of the damage to the cylinder 
1s at the back with some damage to the front and relatively 
little damage to the top and bottom. Reducing the stiffener 
spacing also effects the damage pattern by causing more 
damage to the front of the cylinder and less to the top and 
bottom of the cylinder than for the 12.0 inch sti#rene. 


Sac Ie 
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The resulting damage patterns for a perfect cylinder 
12.0 inches in diameter (shell thickness to radius ratios 
ma@eme 0.0) to 0.04 and stiffener spacing to radius ratios 
meem |.0 to 2.0) are shown in Figures 29 and 30. The shock 
pressure for the 0.06 inch shell thickness is 500 psi for 1 
ms, for the 0.12 inch shell thickness the pressure is 1000 
psi for 1 ms and for the 0.24 inch shell thickness the 
pressure is 1900 psi for 1 ms. 

At a shell thickness of 0.06 inches the pinching of the 
sigetl occurs for stiffener spacing of 12.0. 8.0 and 6.0 
inches. At a stiffener spacing of 12.0 inches increasing 
the shell thickness causes the pinching of the shell to 
decrease until it is eliminated for the shell thickness of 
0.24 inches. At stiffener spacing of less than 12.0 inches 
increasing the shell thickness causes the pinching of the 
shell to disappear for both the 0.12 and 0.24 inch shell 
thickness. 

The end view shown that, for shell thickness of 0.06 
and 0.12 inches, at a diameter of 12.0 inches there is a 
protrusion of the shell toward the shock wave approach 
eemeeeclion. Increasing the shell thickness to 0.24 inches, 
for the three stiffener spacing investigated, causes the 
protrusion to be eliminated. Decreasing the stiffener 


spacing causes a decrease in damage to the cylinder and some 
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minor shifting of the damage pattern but the basic shape of 
the damages pattern remains constant. 

The resulting damage for a perfect cylinder 18.0 inches 
in diameter (shell thickness to radius ratios from 0.0067 to 
0.0267 and stiffener spacing to radius ratios from 0 7G 7c 
1.33) 1s shown in Figures 31 and 32. The shock pressure for 
the 0.06 inch shell thickness is 500 psi for 1 ™sS, forum 
0.12 inch shell thickness the pressure is 1000 psi for 1 ms 
and for the 0.24 inch shell thickness the pressure is 1700 
DSi fers mse 

In all cases, for the 18.0 inch diameter perfect 
cylinder there was no pinching of the shell near the 
stiffener. This cylinder is structurally much less "sium 
than the smaller diameter cylinders thus allowing the shell 
to collapse more readily without pinching near _ the 
stiffener. 

From the end view it can be seen that the protrusion of 
the shell toward the explosion source occurs in all cases. 
For the thicker shell this protrusion is not as sharp asuieem 
the thinner shell. Decreasing the stiffener spacing does 
not greatly affect the damage patterns but does cause the 
damage levels to decrease. 

The damage pattern for a 12.0 inch diameter cylinder, 
with a mode 6 initial imperfection, an amplitude of 5% of 


the shell thickness and no phase shift, is shown in Figures 
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33 and 34. The shock wave pressure varies from 500 psi for 
imiemcor the 0.06 inch shell thickness to 1000 psi for 1 ms 
memeetne O.12 inch shell thickness and finally to 1900 psi 
heme 1 ms for the 0.24 inch shell thickness. These shock 
wave pressure are the same as for the perfect cylinder cases 
at a cylinder diameter of 12.0 inches. 

The front view shows the pinching of the shell near the 
Stiitftener only occurs for the 0.06 inch shell thickness. 
For thicker shells the pinching is eliminated. The end 
views show that the initial mode 6 imperfection 1s seen in 
the final deformation pattern of the shell of the cylinder 
in all cases. Changing the shell thickness and the 
stiffener spacing has very little effect on the damage 
pattern. The damage level is reduced as the stiffener 
Spacing is reduced and this is expected as the cylinder 
becomes stiffer with decreasing stiffener spacing. 

The damage to an 8.0 inch cylinder, with initial 
imperfections consisting of the first 10 mode shapes with an 
amplitude of 1% of the shell thickness and with random phase 
shifts of the modes is shown in Figures 35 and 36. In these 
cases, the shock wave pressures had to be reduced from the 
pressures used in the perfect cylinders at the same 
diameter, because the cylinders completely collapsed at the 
pressures used in the perfect cylinder cases. The shock 


wave pressures were 500 psi for 1ms for the 0.06 inch shell 
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thickness, 1000 psi for 1 ms for the 0.12 inch ema 
thackness and 2250 psi for 1 ms for the 0.24 inch Sopeum 
thickness. At this diameter the cylinder was weakened to a 
greater extent, than for the larger diameter cases, by the 
initial imperfections introduced into the shell geometry. 

The front view of the cylinder shows that there is very 
little pinching of the shell of the cylinder near the 
stiffeners. The end view shows the damage 1s much more 
localized to a smaller portion of the shell than for the 
larger diameter cylinders. The damage pattern is fairly 
consistent with only changes in the magnitude of the damage 
among the cases’. studied. In addition the imam 
imperfection pattern is still seen in the final damage 
pattern but the match is much less pronounced than for the 
later cases studied. 

For a cylinder diameter of 12.0 inches with initial 
imperfections consisting of the first 10 modes shapecui aan 
an amplitude of 1% of the shell thickness with random phase 
shifts of the modes, the damage patterns are much more 
consistent. Figures 37 and 38 shown the resulting damage 
patterns. At this diameter the shock wave pressure is the 
same as for the perfect cylinder cases at the same diameter 
with pressures of 500 psi for 1ms for the 0.06 inch shell 


thickness, 1000 psi for 1 ms for the 0.12 inch @em am 
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wemermmess and 1900 psi for 1 ms for the 0.24 inch shell 
thickness. 

The front view shows a similar damage pattern for all 
cases. The shell very distinctly dishes between the 
Stiffeners in a manner similar to that seen in many 
underwater shock test cylinders. Changing the shell 
thickness and the stiffener spacing has very little effect 
on the damage pattern as seen from the front of the 
cylinder. The end view shows again that in all cases shown 
the damage pattern remains very consistent and only the 
damage levels change with changes in the stiffener spacing. 
The 0.24 inch shell thickness does show some shifting of the 
damage pattern but the same basic pattern is still evident. 

Finally, Figures 39 and 40 show the deformation of an 
18.0 inch diameter cylinder with initial imperfections 
consisting of the first 10 mode shapes, each with an 
amplitude of 1% of the shell thickness with random phase 
Shift. Again, the shock wave pressure was the same as for 
the perfect cylinder cases with a pressure of 500 psi for l 
imewere@r tne 0.06 inch shell thickness, 1000 psi for 1 ms for 
the 0.12 inch shell thickness and 1700 psi for 1 ms for the 
0.24 inch shell thickness. 

From the front view it can be seen that the shell of 
the cylinder has deformed by dishing between the stiffeners 


in all cases. Decreasing the stiffener spacing reduced the 


Sh) 


magnitude of the deformation but did not change the shape of 
the deformation. Increasing the shell thickness had little 
effect on the deformation pattern. The end views of the 
cylinder all show the same basic deformation pattern. 
Increasing the shell thickness had some effect on the damage 
pattern but the same deformation pattern is present in all 


the deformed cylinders. 
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V. THREE DIMENSIONAL RING STIFFENED FINITE LENGTH CYLINDER 


From the three-dimensional infinite cylinder a three 
dimensional finite length cylinder was developed as shown in 
Figure 41. This cylinder is three feet long, 1 foot in 
diameter with two stiffeners evenly spaced 12 inches apart. 
The shell of the cylinder is mild steel 0.06 inches thick. 
The stiffeners are mild steel 0.12 inches thick and 1 inch 
deep. The endplate is HY-100 steel 0.25 inches thick. The 
HY-100 steel is modeled as an elastic-plastic material with 
a Young's modulus of Dy chalie ¢ Wouwmmeteeotsson'’ Ss ratio of 0.3, 
a hardening modulus of 5020 psi and a yield stress of 
moe, 000 psi. This cylinder is modeled as a half cylinder 
with a plane of symmetry perpendicular to the axis of 
rotation. The model has 40 elements in the circumferential 
direction and 15 elements in the axial direction for a total 
of 600 elements and 921 nodes. The use of a half symmetry 
model with a symmetric boundary condition results in a 
smaller number of elements and greater computational 
efficiency for the finite element analysis. For details of 


the model input parameters see Appendix E. 
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Figure 41. Three Dimensional Ring Stiffened Finite Length 
Cylinder Model 


a2 


A. Perfect Cylinder Models 

This cylinder was subjected to the same planar shock 
wave as in the previous models. The resulting damage 
pattern for a perfect cylinder can be seen in Figures 42 
fereugh 49. Again the side of the cylinder facing the 
planar shock wave shown a local raised area that was also 
seen in the previous models. In addition there is a pinch 
in the shell near the endplates and on either side of the 
stiffeners. This pinch was also noted by Chisum (1992) 
during his analysis of DNA models used for underwater shock 
testing. Most of the deformation of the cylinder shell 
occurs between 1.0 and 2.0 milliseconds (Figures 44 through 
ey); PiainigG the tame tsom 1.0°°to 2.0 milliseconds this 
kinetic energy is transferred to strain energy causing 
deformation of the shell. After 2.0 milliseconds the most 
of the plastic deformation has occurred and the cylinder 
shell then deforms only through elastic vibration (Figures 
48 and 49). Figure 50 is a plot of the kinetic energy of 
the cylinder shell as a function of time. The shell has its 
maximum kinetic energy at 1 ms and then this energy rapidly 
decreases as it is dissipated into the surrounding medium 
and into the shell material as strain energy. Figure 51 
Shows the strain history of the element on the front of the 
cylinder exactly in the middle of the stiffeners. From this 


plot, at 1 millisecond the hoop strain is -5590 microstrain, 
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Deformation of Three Dimensional Ring Stiffened 
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Front View @ 1.0 ms (deformations scaled up 2x) 


Figure 44, Deformation of Three Dimensional Ring Stiffened 
Finite Length Cylinder Subjected to a Square Pressure Pulse 
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End View @ 1.0 ms (deformations scaled up 2x) 


Figure AVS) A Deformation of Three Dimensional Ring Stiffened 
Finite Length Cylinder Subjected to a Square Pressure Pulse 
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Front View @ 2.0 ms (deformations scaled up 2X) 


Figure 46. Deformation of Three Dimensional Ring Stiffened 
Finite Length Cylinder Subjected to a Square Pressure Pulse 
of 500 psi for 1 ms 
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End View @ 2.0 ms (deformations scaled up 2X) 
Figure 47. Deformation of Three Dimensional Ring Stiffened 


Finite Length Cylinder Subjected to a Square Pressure Pulse 
Of 500 psi for sl ms 
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Front View @ 5.0 ms (deformations scaled up 2X) 


Figure 48. Deformation of Three Dimensional Ring Stiffened 
Finite Length Cylinder Subjected to a Square Pressure Pulse 
Ge SUC esi. tome.) ms 
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End View @ 5.0 ms (deformations scaled up 2X) 
Figure 49. Deformation of Three Dimensional Ring Stiffened 


Finite Length Cylinder Subjected to a Square Pressure Pulse 
Of 00s ps1 for 1 ms 


> 


600 








S 
=) 400 
oN 
~ 
® 
= 
Lu 
2 
E 200 
Y 
0 
0 .001 002 .003 .004 .005 
Time (sec) 
Figure 50. Kinetic Energy of Shell Material 
0 
£ 
© -0.6x10° 
Y) 
© 
6 
= 
Q. 
© 
eae 
-1.8x104 
0 .001 002 .003 .004 .005 
Time (sec) 


Figure 51. Hoop Strain of Front Middle Element 
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at 2 milliseconds the hoop strain is -15000 microstrain and 
at 5 milliseconds the hoop strain is -14600 microstrain. 
The strain at 5 milliseconds is permanent plastic strain. 
Most of this plastic strain occurred between 1 and 2 
milliseconds with very little additional strain occurring 
after 2 milliseconds. This confirms the previously noted 
timing of the deformation of the shell of the cylinder. 
Most of the damage to the shell occurs between 1 and 2 


Mmelliseconds. 


B. Imperfect Cylinder Models 

The addition of initial imperfections greatly changes 
the resulting deformation pattern of the shell of the 
cylinder. The addition of a mode 6 imperfection (shown in 
Figure 52) with a modal amplitude of 5% of the shell 
thickness and with no random phase shift results in the 
final damage pattern shown in Figure 53 and 54. The end 
view clearly shows the strong effect of the mode 6 initial 
imperfection causing the final deformation to follow this 
initial imperfection. The front wien of the cylinder shows 
that the pinching of the shell near the endplates and 
stiffeners is very much reduced due to the introduction of 
this imperfection. 

An investigation into the effect of the modal amplitude 


on the final deformation pattern showed that very small 
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Initial Imperfection (imperfections scaled up 100X) 


Figure 52. Three Dimensional Ring Stiffened Finite Length 
Cylinder Subjected to a Square Pressure Pulse of 500 psi for 
1 ms, Mode 6 Initial Imperfection As=0.05h no phase shift. 
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Front View @ 5.0 ms 


Figure 53. Three Dimensional Ring Stiffened Finite Length 
Cylinder Subjected to a Square Pressure Pulse of 500 psi for 
1 ms, Mode 6 Initial Imperfection As=0.05h no phase shift 
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End View @ 5.0 ms 
Figure 54. Three Dimensional Ring Stiffened Finite Length 


Cylinder Subjected to a Square Pressure Pulse of 500 psi for 
1 ms, Mode 6 Initial Imperfection As=0.05h4 no phase shift 


61 


modal amplitudes have a significant effect on the response 
of the cylinder. An initial imperfection consisting of the 
first 10 mode shapes with varying modal amplitudes and 
random phase shifts 1s shown in Figure 55. The shell 
deformation again follows the initial imperfections present 
in the shell structure. As the modal amplitudes increase 
the deformation of the shell of the cylinder also increases. 
At a modal amplitude of 0.01% of the shell thickness 
(Figures 56 and 57) the front of the cylinder is diStamemm 
raised toward the explosive charge as in the perfect 
cylinder. But, it can also be seen from the end view that 
the deformation pattern even at this magnitude of 
imperfection has changed from the pattern for the perfect 
cylinder. At a model amplitude of 0.05% (Figures 58 and 59) 
this alteration of the pattern is quite obvious. The 
pinching of the shell of the cylinder persists until higher 
modal amplitudes are reached (Figures 60 through 63). 
Finally at a modal amplitude of 0.5% of the shell thickness 
(Figures 62 and 63) the pinching of the shell is eliminated. 
Increasing the amplitude further results in increased damage 
to the shell due to the shock wave (Figures 64 through 67). 
It should be noted that an amplitude of 0.5% of the shell 
thickness represents a modal amplitude of 0.0003 inches for 
this model. This magnitude of imperfections would almost 


certainly be present in a test cylinder. 
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/ 


Initial Imperfection (imperfections scaled up 200X) 
Figure 55. Three Dimensional Ring Stiffened Finite Length 


Cylinder Subjected to a Square Pressure Pulse of 500 psi for 
1 ms, First 10 Modes with Random Phase Shift 
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Front View @ 5.0 ms An»=0.0001h with Random Phase Shift 

Figure 56. Deformation of Three Dimensional Ring Stiffened 

Finite Length Cylinder Subjected to a Square Pressure Pulse 
of 500 psi for 1 ms 
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End View @ 5.0 ms An=0.0001h with Random Phase Shift 

Figure 57. Deformation of Three Dimensional Ring Stiffened 

Finite Length Cylinder Subjected to a Square Pressure Pulse 
of 500 psi for tems 





Front View @ 5.0 ms An=0.0005h with Random Phase Shift 


Figure 58. Deformation of Three Dimensional Ring Stiffened 
Finite Length Cylinder Subjected to a Square Pressure Pulse 
of 500 psi for 1 ms 
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End View @ 5.0 ms An=0.0005h with Random Phase Shift 
Figure 59. Deformation of Three Dimensional Ring Stiffened 


Finite Length Cylinder Subjected to a Square Pressure Pulse 
GEIS 0C4 Osi, for | ins 
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Front View @ 5.0 ms n=Q0.001k with Random Phase Shift 


Figure 60. Deformation of Three Dimensional Ring Stiffened 
Finite Length Cylinder Subjected to a Square Pressure Pulse 
Of 5U0spSietor) 1 ms 
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End View @ 5.0 ms An=0.001Ih with Random Phase Shift 


Figure 61. Deformation of Three Dimensional Ring Stiffened 
Finite Length Cylinder Subjected to a Square Pressure Pulse 
Of 500 psi tor ims 
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Front View @ 5.0 ms An=0.005h with Random Phase Shift 


Figure 62. Deformation of Three Dimensional Ring Stiffened 
Finite Length Cylinder Subjected to a Square Pressure Pulse 
of 500 psi for 1 ms 
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End View @ 5.0 ms An=0.005h with Random Phase Shift 
Figure 63. Deformation of Three Dimensional Ring Stiffened 


Finite Length Cylinder Subjected to a Square Pressure Pulse 
Gin 500 Dez tor. 1 ms 
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Front View @ 5.0 ms An»a=0.0Ih with Random Phase Shift 


Figure 64. Deformation of Three Dimensional Ring Stiffened 
Finite Length Cylinder Subjected to a Square Pressure Pulse 


Of 500 psi, tor Lams 
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End View @ 5.0 ms An»aA=0.0lIh with Random Phase Shift 


Figure 65. Deformation of Three Dimensional Ring Stiffened 
Finite Length Cylinder Subjected to a Square Pressure Pulse 


Of 500 Osi vEor ems 
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Front View @ 5.0 ms n=0.05h with Random Phase Shift 


Figure 66. Deformation of Three Dimensional Ring Stiffened 
Finite Length Cylinder Subjected to a Square Pressure Pulse 
of 500 psi for 1 ms 
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End View @ 5.0 ms n=Q.05h with Random Phase Shift 
Figure 67. Deformation of Three Dimensional Ring Stiffened 


Finite Length Cylinder Subjected to a Square Pressure Pulse 
Spel 00) psi for if ms 
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The development of a raised section on the cylinder 
facing the shock wave was unexpected. With the previous 
perfect cylinder the shock wave impacted the cylinder along 
the line of nodes at the front of the cylinder. The 
cylinder was then rotated so that the shock wave would 
impact the cylinder along a line of elements. The resulting 
deformation pattern is shown in Figure 68. There is still a 
local raised section on the cylinder facing the shock wave 
direction. Thus, it is felt that the shot geometry of the 
finite element mesh is not a contributing factor to the 
development of this raised section. 

ge re¥elehuies ejal- a very specific imperfection was 
introduced in the five nodes facing the explosive charge 
centrally located in the midbay of the cylinder. This 
imperfection was an inward imperfection of the node 
locations with an amplitude of 1% of the shell thickness 
(0.006 inches). The damage pattern resulting from a 500 psi 
1 millisecond shock wave is shown in Figure 69. 

The use of Hughes-Liu shell element (Hughes 1981) vice 
the Belytschko-Tsay shell element was investigated. The 
Hughes-Liu shell element does not use many of the 
Simplifying assumptions used in formulating the Belytschko- 
Tsay shell element. As a result, the use of the Hughes-Liu 
Shell element requires longer computational times that may 


not be acceptable in certain models. Figures 70 and 71 show 
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Shock 
Wave 





End View @ 5.0 ms 


Figure 68. Deformation of Three Dimensional Ring Stiffened 
Finite Length Cylinder Subjected to a Square Pressure Pulse 
of 500 psi for 1ms, Cylinder Rotated 4.5 Degrees 





li 
al 







I 


Xx 
\« 
i 
St 
Hy 
a 
Of, 
Oy 


SS 
we 
[~ 
Se, 
Rt 
aS 
Rae 
Hy 
gor. 
HE 
S 
aS 
= 
a 


= 
+ 
Y 
N 
a 
ae 
i 
+H 
Nuuasy 
TH 
ae 
es 
a 
= 
-—| 


Shock 
Wave 


oe 

wees 
HS 
ue 

me 

an 
ae 










SSS Ss 
SSeS ae 


End View @ 5.0 ms 
Begure 69. Three Dimensional Finite Cylinder Subjected to a 


Square Pressure Pulse of 500 psi 1ms, Initial Imperfection 
Front Midbay Nodes 1% Offset Radially Inward 
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the resulting deformation to a model using Hughes-Liu shell 
elements. This model is identical to the previous models 
with the exception of the use of the Hughes-Liu_ shell 
element. This model has initial imperfections with the 
first 10 modes and modal amplitudes of 1% of the shell 
thickness and random phase shifts. Comparison of Figures 70 
and 71 with Figures 64 and 65 shows that the Hughes-Liu and 
the Belytschko-Tsay shell elements have the same final 
deformations. In this case the use of the numerically more 


complicated Hughes-Liu shell element is not warranted. 


C. Exponential Decay Shock Wave 

The use of an explosive shock wave modeled as a square 
wave 1S representative of a very large explosion at a great 
standoff distance from the cylinder. In order to model a 
smaller explosion close to the cylinder a pressure profile 
corresponding to that produced by 40 pounds of PETN at a 
standoff distance of 30 feet was used. The pressure time 
history of the resulting shock wave is shown in Figure 72. 
This explosion produces a peak pressure of 1839 psi and an 
exponentially decaying pressure history with a time constant 
of 0.311 milliseconds. The resulting damage to the cylinder 
from this shock wave is shown in Figure 73 through 75. From 
the offset view it can be seen that the response of the 


cylinder is different from that to the previous plane wave. 
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Front View @ 5.0 ms 


Figure 70. Three Dimensional Ring Stiffened Finite Length 
Cylinder Hughes-Liu Shell Elements Subjected to a Square 
Peessure Pulse of 500 psi for 1 ms, Initial Imperfection 
First 10 Modes An»=0.0lh with Random Phase Shift 
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End View @ 5.0 ms 


Figure 71. Three Dimensional Ring Stiffened Finite Length 
Cylinder Hughes-Liu Shell Elements Subjected to a Square 
Pressure Pulse of 500 psi for 1ms, Initial Imperfection 
First 10 Modes An=0.0lh with Random Phase Shift 
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Shock Wave Pressure (psi) 
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Figure 72. Pressure Time History of Shock Wave from 40 lbs 


PETN with a 30 foot Standoff and 1ms Surface Cutoff Time 
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Front View @ 5.0 ms 
Figure 73. Three Dimensional Ring Stiffened Finite Length 


Cylinder Subjected to a Pressure Pulse from 40 lbs PETN 30 
FOOE Standerse 
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Hagure 74. 


Shock 
Wave 


Figure 75. 





Offset View @ 5.0 ms 


Three Dimensional Ring Stiffened Finite Length 
Cylinder Subjected to a Pressure Pulse from 40 lbs PETN 30 


Foot Standoff 
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End View @ 5.0 ms 


Three Dimensional Ring Stiffened Finite Length 


Cylinder Subjected to a Pressure Pulse from 40 lbs PETN 30 


Foot Standoff Final 
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The small standoff distance causes the shock wave to be a 
spherical wave at the cylinder. The midbay still has a 
protrusion of the shell material toward the explosive charge 
however the endbay shell material has a depression of the 
Shell facing the charge. Also evident from the front view 
1s the severe pinching of the shell near the endplate and on 
either side of the stiffener. This pinching was very severe 
near the endplate. 

The introduction of a 5% mode 6 imperfection to the 
model (Figure 76) results in the deformations shown in 
Figure 77 and 78. Again, it 1s evident that the damage 
pattern of the cylinder due to the explosive shock clearly 
followed the initial imperfection. The pinching GEaieee 
shell material near the endplate and stiffeners was reduced 
but not eliminated. 

Introducing an imperfection of the first 10 mode shapes 
with a modal amplitude of 1% of the shell thickness with 
random phase shifts (Figure 79) results in the final damage 
pattern shown in Figures 80 and 81. The deformation of the 
Shell of the cylinder was again Seneca by the initzal 
imperfection pattern. The local pinching of the shell 
material on either side of the stiffener was eliminated and 


the pinching near the endplate was reduced in magnitude. 
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Initial Imperfection (imperfections scaled up 100X) 


Figure Wise auinaaee Dimensional Ring Stiffened Finite Length 

Cylinder Subjected to a Pressure Pulse from 40 lbs PETN 30 

Foot Standoff, Initial Mode 6 Imperfection As=0.05h no phase 
shift 
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Three Dimensional Ring Stiffened Finite Length 


Figure 77. 
Cylinder Subjected to a Pressure Pulse from 40 lbs PETN 30 
Initial Mode 6 Imperfection As=0.05h no phase 
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End View @ 5.0 ms 


Three Dimensional Ring Stiffened Finite Length 


Figure 78. 
Cylinder Subjected to a Pressure Pulse from 40 lbs PETN 30 
Initial Mode 6 Imperfection As=0.05h no phase 


Foot Standoff, 
Slag 


78 


Initial Imperfection (imperfections scaled up 200X) 


Figure 79. Three Dimensional Ring Stiffened Finite Length 

Cylinder Subjected to a Pressure Pulse from 40 lbs PETN 30 

Foot Standoff, Initial Imperfection First 10 modes A,»=0.0Ih 
with Random Phase Shift 


US 






| 


ng 












nae oo eee 
See 





it 
\ 
1 


Front View @ 5.0 ms 


Figure 80. Three Dimensional Ring Stiffened Finite Length 

Cylinder Subjected to a Pressure Pulse from 40 lbs PETN 30 

Foot Standoff, Initial Imperfection First 10 modes A.»=0.01h 
with Random Phase Shift 
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End View @ 5.0 ms 


Figure 81. Three Dimensional Ring Stiffened Finite Length 

Cylinder Subjected to a Pressure Pulse from 40 lbs PETN 30 

Foot Standoff, Initial Imperfection First 10 modes An=0.01h 
with Random Phase Shift 


80 


VI. SUMMARY AND CONCLUSIONS 


The response of test cylinders to underwater explosive 
shock 1S a complicated function of many factors. One 
important factor that is often overlooked in the modeling of 
these underwater shock phenomena is the Dds ee 
imperfections that are present in the test cylinders. The 
response of model cylinders Subjected to Simulated 
underwater explosive shock has been shown to be very 
dependent to the initial imperfections introduced into the 
cylinder models. The introduction of these imperfections 
not only caused the shape of the shell of the cylinder (as 
viewed from the end) to follow the shape of the initial 
imperfection but it has also changed the response of the 
Shell near the endplates and stiffeners of the model 
ev lander . The pinching of the shell near these stiffeners 
and endplates was greatly reduced or eliminated by the 
introduction of initial imperfections in the model 
cylinders. This resulting response is closer to the 
response observed in test cylinders subjected to actual 
underwater shock loading. im ocdition, , Ehe initial 
imperfections introduced increased the magnitude of the 
Shell deformations compared to the deformations’ for 


cylinders modeled as perfect cylinders. 
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Changes in the cylinder geometry resulted in some 
changes in the response of the cylinder to underwater shock. 
However in all cases the InMewecicrk ten.» OF initial 
imperfections greatly affected the response of the cylinder 
as compared to the response of a perfect cylinder. 

If the initial imperfections of a test cylinder are 
known and introduced into a model cylinder used in finite 
element analysis, the results of the finite element analysis 
may more closely, simulate the actual response of the test 


cylinder. 





APPENDIX A 


FORTRAN PROGRAMS FOR MODIFYING INGRIDO FILE FOR MODAL 
TMEER Eee tlTONs 


program imp 
KKK KK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK 
lees, OrOOram written by Donald T. Hooker II on 3/15/93 * 
Revised 5/12/93 e 
KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEKE KKK KEKE KKK KEK 
This program modifies the nodes of a file called 
"“ingrido.raw", which is an edited version of ingrido in 
which the first line is the beginning of the node list. 
This 1S one of a series of 2 programs, which need to be 
run in order, 1.e. "imp" then "convert". These programs 
modify the node positions for modal imperfections WITH 
OR WITHOUT random amplitudes and/or phase shifts. The 
final output file name is imp.dat. 
KR KK KKK KKK KEKE KKK KKK KEK KK KKK KK KKK KK KKK KKK KK KKK KEKE KKK KKK KEKE KKK KK 
VARIABLE DECLARATION 
a - angle from reference axis to radial position of 
node 
al thru alQ - weighting factor for each modal 
imperfection 
bl - ingrido constants 
b2 - ingrido constants 
dr - change in radial position due to modal 
imperfections 
drl thru drl10 - radial position change for each modal 
imperfection 
h - thickness of shell 
n - node number 
m - number of nodes 
eee 35 ..14159265359 
psl thru ps10 - phase shift for each modal 
imperfection 
r - radial position of node 
randamp - if random amplitude is desired set randamp=1* 
randphs - if random phase is desired set randphs=1 
hn x COOLGdinate of node 
Y - y coordinate of node 
Z- z coordinate of node 


KKK KKK KKK KEKE KKK KKK KKK KKK KKK KEKE KKK KKK KKK KKK KEKE KEKE KEK KKK KEKE KKK KEKE 


* 


+ + + £ £ + £ £ FF € € FF £ FF FF £ iF F: «iF OF FF: F OF * OF ~ OF ~* OF 


* 


* 


* 


+ + + + + ££ + + FF ££ ££ FF FF FF FF iF F OF ~*~ F ~*~ iF ~*~ OF ~* OF ~*~ OF:~* OF * OF :~ OF KH OF ~ OF ~ iF ~« O * K€ OK OF 


G@OWewe Precision Xy7,2Z,a,r)ar;h;p1 
double precision al,a2,a3,a4,a5,a6,a7,a8,a9,al10 
Geotmeme Precicton Gril, dr2,dr3,dr4,dr5,dr6,dr/,drée,dar9, 
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aoa di), 

double precision psl,ps2,ps3,ps4,ps5,ps6,ps7,ps8,ps9, 
& psl0 

integer randamp, randphs 


KEK KEKE KEKE KK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KK KEK KKK KKK KK 


is OPEN INPUT AND OUUEUT Fires es 


KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKK KEK KKK 


epen (14,7 file="1nemide naw 
open(15, file='imp.unconverted' ) 


KKKKK KKK KKK KKK KEK KKK KKK KKK KKK KK KKK KK KKK KK KKK KKK KKK KK KKK KK KKK 


Li IF RANDOM MODE AMPLITUDE IS DESIRED SET RANDAMP EQUAL* 
= TO. % 
J IF RANDOM MODE PHASE SHIFT IS DESIRED SEI Ante. i 
is EQUAL TO 1 ig 


KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KK KKK KKK KKK KKK KKK KKK KKK 


hamcaneoe=—— 0 

randphs = 1 
KKK KKK KEKE KKK KKK KKK KKK KKK KKK KEK KKK KK KKK KK KKK KKK KKK KKK KKK KKK KKK 
- INPUT NUMBER OF NODES TO BE MODIFIED AND WRITE NUMESeS 
is TOVOUTEUT FILE Tie UNC eh ERageie 


KKK KKK KKK KEKE KKK KKK KKK KKK KKK KKK KKK KKK KEK KKK KK KKK KKK KKK KKK 


Wee Seal 
write(15,98)m 
98 format (<1 5) 


KKK KKK KEKE KKK KEKE KK KKK KKK KKK KKK KKK KKK K KKK KKK KKK KK KKK KK KKK KKK 


INPUT WR TGKNES Seer. on eEe = 


KKK KEKE KKK KEKE KKK KK KEK KKK KKK KEKE KEKE KEKE KK KKK KKK KKK KKK KKK KK KKK KKK KK 


levee 0) 


KEK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKK K KKK KK KKKKKKK KKK KKK 


. SBT THE VALUE VOR pri 


KKKKKK KKK KK KKK KKK KEK KEKE KKK KEKE KKK KEKKK KKK KK KKK KKK KK KKK KKK KKK KKKK 


Ol 3145 9265 59 


KR KKK KEKE KKK KKK KKK KK KK KKK KKK KKK KK KKK KKK KKK KKK KKK KK KKK KKK KKK KKK 


ig INPUT THE WEIGHTING COEFPICUTENTS VOR Wy rete couiaman | is 
a MODE SHAPES ig 


KKKKKKKKKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KK KKK KKK KKK KKK KKK KKK 


a = OO 
AZ = 710 Ou 
sib =O) 2 ON Tl 
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a4] 0201 
aioe — G20 
eor = O01 
av = 0.49 
moe 0) 64 
a7 = 0.81 
eng. = 120 


KKK KEKE KKK KK KEKE KKK KEK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEK KKK KKK KKK KKK 


% MVE PHASE SHIFT FOR BACH MODE SHAPE a 


KEKE KK KEK KEKE KKK KEK KKK KEK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK 


mel = 0.0 
mez = 0.0 
jos = 0.0 
ps4 = 0.0 
poo = 0.0 
mso = 0:0 
es, = 0.0 
eos = 0.0 
poo = 0.0 
slo = 0.0 


KEKE KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KEKE KK KEKE KK KEKKK KKK KKK 


e INPUT SEED NUMBER FOR RANDOM NUMBER GENERATOR - 


MI can KKK Kk eK KK Keke Kae KKK KKK KK KKK KEKE KKK KKK KEKE KKK KK KKK KKK 


yeed = 1.0 
KEK KKK KEKE KKK KEK KKK KEK KEKE KK KEK KEKE KKK KKK KK KKK KKK KEKE KEKE KKK KKK KKK KKK 
id CALL SUBROUTINE SRAND TO GENERATE SEED NUMBER FOR 5 
¢ RANDOM NUMBER GENERATOR SUBROUTINE i 


KEK K KKK KKK KR KKK KKK KKK KKK KKK KKK KE KKK KK KKK KK KKK KKK KKK KKK KKK KKK 


call srand(seed) 


KEK KK KEK KEK KKK KKK KKK KKK KKK KKK KKK KKK KK KEK KKK KKK KEKE KEKKKEKEK KEKE KK 


e Sab SUBROUTINE RAND TO GENERATE RANDOM PHASE SHIFT * 
‘s FOR EACH MODE SHAPE ts 


KEKE KKK KEK KEKE KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEKK KK KKK KK KKK KK 


1f (randphs .eq.1) then 


Poe 2 Boe rand | ) 
sc =.2.. % pa rand) 
eseme2s * DL * rand (} 
ps4 Zee 2 “eamel) 
Geo aez. © pa -* rand () 
eae ete On * rane |) 
ese 224" pr, * rand) 
ecw =—= io. @o2- * rand() 
Bevo we wots” rand) 
Boca Or. * rend () 


endif 


KKK KKK KK KKK KEKE KEK KKK KEKE KEK KKK KKK KK KKK KKK KEKE KE KEKE KEK KK KKK KKK KKK KKK KKK 


i CALL SUBROUTINE RAND TO GENERATE RANDOM WEIGHTING OF * 
i MODE SHAPES : 


KEKKKKKEKEKKKEKKEKREKKEKEKEKEKEEKKREKEKKEKKEKEKKKKKKEKKEKEK KEKE RE ee 


1f (randamp .eq. 1) then 
al = al * rand() 
a2 = a2 * rand() 
a3: =7as  * wane 
a4 = a4 * wand) 
a> =fa54 = surance 
a6 = a6 * rand() 
a?’ = af)? strane 
ag<= ag ~ rand) 
a9 = a9 * rand() 
alQc=.4a b0) * “rand 

endif 


KKK KEK KEKE KEKE KEK KEK KK KEKE KKK KEKE KKK KEK KEKE KEKE KK KEKE KEKE KEKE KEK KEKE KEKE KEKE KK KEKE KEKE KEK 


ig DO LOOP TO READ INPUT DATA ic 


KKK KEKE KEKE KEKE KEKE KEK KKK KEKE KEKE KKK KEK KKK KKK KEKE KKK KEKE KKK KEKE KEKE KEKE KEKE KKK KEKE KKK 


Cle? IO at Sign 
réad( 145) sar ole yar 


KKK KKK KEKE KKK KKK KEKE KKK KEKE KKK KK KEKE KKK KEKE KKK KKK KEKE KEKE KK KEKE KEKE KKK KKK KKK 


ia CALCULATE (RADIUS OF FNODE  FOSl Tren ig 
(Look carefully at the axis of rotation of the i 
. cylinder to determine the proper Cartesian * 
* coordinates used to calculate the radius of the bs 
bs cylinder. In this case the axis of rotation 1S emewz. 
- axis. = 
KKK KEK KKK KKK KEK KKK KKK KKK KKK KKK KKK KKK KKK KK KEK KK KK KEKE KKK KKK KK KKK 


Pa (Ke Sys es 


KEK KEKE KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEKE KKK KEK KKK KKK 
i CALCULATE THE ANGLE FROM THE REFERENCE AXIS TO THE a 
* NODE POINTS (Look carefully at the zero angle ig 
is direction. In this case the x axis points in the . 
: direction of the charge.) a 
KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEK KKK KKK KKK KKK KKK KKK KKK KKK KKK 


1£ ( (dabs(y) Ue. 07001) Seen 
Ges oecohe SO 0) wlayeta 
a= OO 
else 
a=p1 
endif 
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else 
a=datan(y/x) 
timex tte 0.0) then 


a=at+pl 
endif 
endif 


KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEKE KKK KEKE KEKE KEKE KKK KK KKK KEKE 


= Mebiry THE NODE POSITIONS USING THE FIRST TEN MODE S 
ig SHAPES : 


KKK KKK KEK KEK KKK KEKE KK KK KKK KEK KEK KEK KKK KKK KKK KEKE KEK KEKE KK KKK KKKKEKEKEKE KEE 


diala= aly * in * dees la +) ps!) 

ena, = sar Sn = deost( 2.7 a + ps2) 

Giese = ce = ie telcos (3. * a + ps3) 

Giese = ae i= deost4. * a + psd) 

Guam =. a5." ) *sdeos(s. *. a + ps5) 

@eon="ao-* Hh =“ Veeosi6. * a + ps6) 

@iei = a7. * Wf 727/-7.° * deos(/. * 4a + ps7) 

aigoe = tes" fey 75. 7 6. (* -deos (Ss * Ta. ps5) 

elgg = ao = he/ Yue 9. * dcos(9.- * a + psy) 

@iateo) = vol. hn 7 Oe of “10. * Adcos(i0.* a + psl0) 
Glam ete lee ce? Prods -— ard = dro co dre +4) 0r/ =+ “drs 4 


& ar9 + drl0 


KKK KKK KKK KK KKK KKK KKK KKK KKK KE KKK KEKE KEK KKK KKK KKK KK KEKE KEKE KEKE 


J CALCULATE THE CARTESIAN COORDINATES OF THE NODES ie 


KKK KK KKK KKK KR KKK KK KKK KR KEK KKK KKK KKK KKK KKK KEK KEKE KKK KKK KE KKK KEKE KKK KKEE 


x=x-dr*dcos (a) 
y=y-dr*dsin(a) 


KEK KK KEK KEK KKK KKK KKK KEK KKK KKK KKK KKK KEKE KKK KEKE KK KEKEKEKEKEKEKEKKKEKEKEKEKKE KK 


is WRITE DATA TO TEMPORARY DATA FILE IMP.UNCONVERTED e 


KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KEKE KEKKE KEKE KE KE 


Wiatee( 5,998) n,bl,x,Vv,Z, 2 
998 meamat(1s,f5,0,3e20.13,£5.0) 
100 continue 


KKK KKK KK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEKE KEKEKKEKKKKEKEEK 


Fp @peoe INPUY AND OUTPUT FILES & 


KKK KK KKK KK KKK KKK KKK KK KKK KKK KK KKK KEK KEK KEKE KEKE KKK KKK KEKE KKK KKEEKEEE 


close(14) 
close(15) 
stop 

end 
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program convertimp 


KKK KK KR KKK KKK KKK KR KRAKKKAKK KKK KKKKKEKEKRERKKKRARKKKEKKEKR KEKE Ke 


* This program written by Donald T. Hooker Ii on 3/15/7230 
* Revised 3/30/93 . 
KKK KKK KKK KKK KKK KEK KKK RK KKK KKK KK RK KK KK KKK KKK KKK KKK KKK KKK KK KK 
* This program converts numerical data 
from"imp.unconverted"* 

* anto ascii format for use in modifying ingrido for me@deres 
* imperfections re 


KEKE KKK KEK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEKE KEKE KKK KKK KKK K KKK 


* VARIABLE DECLARATION iy 
* f* - data to be converted to ascii format bs 
* m - number of nodes . 


kaxekkkkkekKekekeKe KKK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KR KKK KKKKRKKKKKEEKEE 


Character *i2 aim 
character* i] f24r3, iz 
character*18 £4 
Charaéwer* “£57 7£6 
character*18 £7 
Characeer* 1 siaestg 
Character*18 £10 
Character sab ial 
open(14, file='imp.unconverted' ) 
open (15 ,flle= oameecdats) 
ff Ze 


kKaeK KKK KKK KKK KK KKK KR KKK KK KKK KKK KKK KKK KKK KK KR KKK KKK KKKRKKK KKK KEK 


S READ NUMBER OF NODES CeeE MeCbIeTes % 


KKK KKK KKK KEK KEK KKK KKK KKK KKK KKK KKK KK KEK KKK KEK KKK KKK KEK KKK KKK KK KK 


read (1477) m 


KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEK KKK KKK KEK KKK K KK KKK KKK KKK 


iy DO LOOP TO READ DATA AND CONVERT IT TO ASCII FORMA. 


KKK KKK KKK KKK KKK KK KEK KKK KKK KEKE KKK KEK KKK KEKE KKK KKK KKK KKK KKK KKK KKK KK 


Con 0a— ie im 
read(14,998)f£1,£2,£3,£4,£5, £67077 20, 27, 2a ae 
998 format (al3,al,al,al8,al,al,al8,al,al,al8,a5) 
write(15,998)£1,£3,£2,£4, 16 7t262) £9, £2 2 Oe 
100 continue 


kKkkkKk ka K Ke KK RK KKK KKK KKK KKEKRKKKKKRKRKKKKKR KK KKKKKRKRKRKKKKKKEKKEE EAS 


CLOSE INPUT AND OUTPUT FILES : 
KEKE KKK KKK KKK KKK KK KKK KKK KK KK KKK KKK KKK KKK KK KK KKK KKK KK KKK kK 
close(14) 
close( 15) 
SEOp 
end 


88 


APPENDIX B 


PicstD INPUT FILE FOR TWO DIMENSIONAL INFINITE CYLINDER 
MODEL 


Piette CYLINDER MODEL 
dn3d vec term 5.0e-3 plti 1.0e-5 prti 1000. 


mat 1 type 3 e 2.9e+7 pr 0.3 ro 7.356e-4 
etan 5.1e+4 sigy 3.2e+4 shell quad 5 thick 0.060 endmat 


fess 0.0 0.0 1.0 0.0 
ieee 2 0.0 2.5e-8 1.0 2.5e-8 
Elan 2 
wo, 003 O 0-1 0.001 symm 
oO 003 Oo loool sya 
Start 
~l1 6 -1ll1; 
-] 6 -1l1 ? 
i - 
-l. on i 
-1. On ih 
=7.00 3 mos 
ea@ewerd 3 3 0 3 6.0 
emi —-S - -1 -3 >; ; 1 -1.0 0. O. O. 
mate 1 
end 
end 
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APPENDIX C 


MATERIAL MODEL DETAILS 


For the steels used in the numerical models, a bilinear 
Stress-strain type model was used. The material model used 
was a "Kinematic/Isotropic Elastic Plastic" model (Hallquist 
1990). Figure 82 shows a typical stress-strain diagram for 
the materials used. A hardening parameter of 0, giving pure 
kinematic hardening, was used for all of the steels. Strain 
rate hardening was not used for the materials in this study. 

Typical properties of the materials used in the 


numerical! models are: 


Mild Steel HY-100 
Density (p) (lb/ft?) : 490 490 
Poission's Ratio(v) Oe O28 
Yield Stress (oy) (psi) : 3.2104 1.08x10° 
Young's Modulus (E) (psi) : 2.9x10/ 2.9x10 
Hardening Modulus (Eran) (psi) ;: 5.1x10° 5. O2sO= 


For the deformations encountered in this study the materials 


did not approach their fracture stress. 
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Stress (ksi) 


125 
100 
ths, 
90 


25 tt 


Prise C2. 





—-— HY-100 Steel 
Mild Stee! 


Strain 


Typical Bilinear Stress-Strain Diagram 
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APPENDIX D 
INGRID INPUT FILE FOR THREE DIMENSIONAL STIFFENED INET 
Cyr Eee Mera 


STIFFENED INFINGD TE CY bINBER MesEE 
dn3d vec term 5.0e-3 plti 1.0e-5 prti 1000. 


Material definitions 


shell material (mild steel) 


O-aee 0 


mat 1 type 3 e 2.9e+7 pr 0.3 ro 7.356e-4 
etan 5.1e+4 sigy 3.2e+4 shell quad 5 thick 0.060 endmat 


& stiffener material (mild steel) 
e 


mat 2 type 3 e 2.9e+7 pr 0.3 ro 7.356e-4 
etan 5.1e+4 sigy 3.2e+4 shell quad 5 thick 0.12 endmat 


Cc 
ec load curve definitions 
e 
led Av -0 S070 0 i Ome 
led 2°2 020 12-06 =6 170 1.0e-6 
S 
e; symmetry planes 
c 
plan 2 
0 0 -6. 0 0 =-1 702001) sy am 
O20 ase 0 0 LF C200 Ssvim 
G 
@ construct shell 
e: 
Stalre 
-1 11 -21; 
=]! = 2a 
lL. otG Hee: 
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= || 0 ae 

ieee O.)6|OCOUl. 

aoe O. 6. 

pan. 3 3). 0 C3. 6.0 

frame) —-3 + -1 -3 > ; 1 -1.0 O. 
mate 1 


end 


G 
S 
(Ss 


S 


il 
IL 


Somstruct Stiffener 
Cart 


pee 22 23 ; 
ome 22 23; 
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OF 


ue 


APPENDIX E 


INGRID INPUT FILE FOR THREE DIMENSIONAL STIFFENED F tae 
CYLINDER MOLEL 


STIFFENED 7 ENTE CY ENE 
dn3aq vec term 5.0e-3 plta” 1406-5 ori Loo 


G material definitions 
e shell material (mild steel) 


mat 1 type 3 e 2.9e+7 pr 0.3 ro 7.356e-4 
etan 5.1le+4 sigy 3.2e+4 shell quad 5 thick 0.060 endmat 
e 
C stiffener material (mild steel) 
@ 


mat 2 type 3 e@ 2.9e47 pr 043 %eo 7.356e-4 
etan 5.1e+4 sigy 3.2e+4 shell quad 5 thick 0.12 endmat 


@ endplate material (HY-100 steel) 
Mat..3cbype 3.6 2e9er) ore). jar oe, olen 

etan 5.02e+4 sigy 1.08e+5 shell quad 5 thick 0.25 
endmat 


e 
é 16ad Curve: Get inue ons 


GC symmetry planes 
Cc 


pian 1 
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a Gn; 
-] on 6 
-1 0.) <1.. 
oO” 0 she 
eee OO 5 35 0 3 6.0 
Meee -l -3 ; -1 -3 >; ; 1 =-1.0 O. 
mate 1 
end 
e 
c Senrseruct stiffener 
(@ 
start 
ieee? 12 13 ; 
fees el 2 13 CS; 
-l ; 
-1 -1 011 
-l1 -1 0141 
6.0 
emus 0 45 ;12045; ; 
cme = l6hlUGShlUCOOtC SC LO 
oo 10) A 6 RR io 
opens 0 4 4 OQ 
mate 2 
end 
e 
e construct endplate 
e surface definitions 
e 
fememevta 0 00 001 6.0 
Somers 0 00 001 ([6.0*3/5] 


fewer 15 19 ; 
ieee 15 19 ; 


ee 
wie) 1 1 
oe OO 1 1 
bey ans 


OF 


Oe 


mee 5 UhlU El LC hULdtlCLCUK1L LO CUO 0.0 0.0 


Pee 4 bls cdl vc 204 5G 


Siva -2 -4 > -2 -4 : +; sd 2 
sfi -1 -5 ; -1 -5 ; ; sdl 
mate 3 

end 

end 


Je 
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